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ABSTRACT 

The Kepler mission has yielded a large number of planet candidates from among the Kepler Objects 
of Interest (KOIs), but spectroscopic follo-w-up of these relatively faint stars is a serious bottleneck in 
confirming and characterizing these systems. We present motivation and survey design for an ongoing 
project with the SDSS-HI multiplexed APOGEE near-infrared spectrograph to monitor hundreds of 
KOI host stars. We report some of our first results using representative targets from our sample, which 
include current planet candidates that we find to be false positives, as well as candidates listed as false 
positives that we do not find to be spectroscopic binaries. With this survey, KOI hosts are observed 
over ^ 20 epochs at a radial velocity precision of 100 — 200 ms“'. These observations can easily 
identify a majority of false positives caused by physically-associated stellar or substellar binaries, and 
in many cases, fully characterize their orbits. We demonstrate that APOGEE is capable of achieving 
RV precision at the 100 — 200 m s“' level over long time baselines, and that APOGEE’S multiplexing 
capability makes it substantially more efficient at identifying false positives due to binaries than 
other single-object spectrographs working to confirm KOIs as planets. These APOGEE RVs enable 
ancillary science projects, such as studies of fundamental stellar astrophysics or intrinsically rare 
substellar companions. The coadded APOGEE spectra can be used to derive stellar properties (Teu, 
log 5 ) and chemical abundances of over a dozen elements to probe correlations of planet properties 
with individual elemental abundances. 
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1. INTRODUCTION 
1 .1. Kepler’s Planet Candidates 

The Kepler spacecraft’s primary mission is to deter¬ 
mine the frequency of Earth-sized exopla nets orbiting in 
the habitable zone of the ir parent stars (iBorucki et al. I 
120101: iKoch et al. I l2010t) . with a second objective of 
studying a wide variet y of stellar ast r ophys ics via as- 
teroseismology (e.g., IChaplin et al. ] [200 . In ad¬ 
dition, the high precision photometry 80 ppm 
over 6-ho ur timescales for the brightest (Kr, S. 15) 
dwarfs, [CaldweH eEaLj 120101 : iGilliland et al. I [201 it 
iChristiansen et al. Il20ldi enables studies of giant exo¬ 
planets and a wide variety of variable stars. Its photo¬ 
metric band Kp covers 423 - 897 nm and is si milar to, but 
broad er than, a combined V and R band (|Koch et al. I 
l2010f) . To find exoplanets, Kepler makes use of the tran¬ 
sit method, which detects planet candidates by mea¬ 
suring the flux loss that occurs when a planet crosses 
the face of its parent star. However, there are several 
sources of false positives that must be taken into account 
when analyzing these candidates, most notably: graz¬ 
ing eclipsing binaries (EBs), EBs (including hierarchical 
triples) whose eclipse depths are diluted by another star 
through flux contamination, brown dwarfs or low mass 
stars that have radii comparable to giant exoplanets, and 
even larger exoplanets that transit a fainter star within 
the photometric aperture. 

Because of these sources of false positives, the Kepler 
team makes a very clear distinction between candidate 
exoplanets and those that have been dynamically con¬ 
firmed through spectroscopic radial velocity (RV) mea- 
surements or through photodynamical m odeling (e.g., 
iHolman et al. I I 2 OIOI : [Carter et aL~l 1201111 . Kepler Ob¬ 
jects of Interest (KOIs) consist of candidate exoplanets, 
eclipsing binaries, and known false positives. Those KOIs 
that are not known to be false positives or EBs are re¬ 
ferred to as “active planet can didates” (jBorncki et al. I 
l2011alibl : iBatalha et al. I l2013il . but for simplicity, we 
will refer to such Kepler planet candidates as “KPCs” 
throughout the rest of this paper. An intermediate level 
of classification consists of “validated” exoplanets, which 
have very low probabilities of being blended EBs as de¬ 
termined through a Monte C arlo statistical anal ysis of 
the Kepler photometry (e.g., [Torres et al. ll2f)TTh . 

As of October 2014, there are a total of 4229 KPCs 
amongst 3251 Kepler stars Q , but only 20% (653) of 
the stars host multiple KPCs. It is estimated that as 
many as 15-26% of transiting plane ts may have clearl y 
detected transit timing variations (jEord et al. 1 1201211 . 
which allow for mass determinations photometrically. 
Even still, a majority of KPCs will require RV obser¬ 
vations to confirm their planetary nature. Such time- 
series RV observations are resource intensive, so effi¬ 
cient identification of false positive candidates is neces¬ 
sary to ensure efficient follow-up of likely planets. In 
addition to aiding in the confirmation of KPCs, robustly 
determining the false positive rate amongst KPCs is re¬ 
quired when conducting statistical analyses of this pop¬ 
ulation. A number of studies have attempted to per¬ 
form such analyses, including investigations of planet 

^ http://exoplanetarchive.ipac.caltech.edu/cgi-bin/ 

ExoTables/nph-exotbls?dataset=ctmiulative_oiily 


frequency as fu nctions of orbital peri o ds and stellar 
host properties (IBorucki et al. I l2011bt lYoudin I 120111 : 
[Howard et al. Il2012il. and studies of the eccentricity dis- 
tribution ([Moorhead et al. Il201lll . 

Aside from the false positive rate of KPCs, knowledge 
of the host star(s) intrinsic properties (e.g., mass, radius, 
effective temperature, surface gravity, metallicity) is re¬ 
quired to determine the masses and radii of the exoplan¬ 
ets, as well as to conduct studies of planetary properties 
as functions of thes e stellar parameters . The Kepler In¬ 
put Catalog (KIC, iBrown et al. ll20lTi l provides a pho¬ 
tometrically derived Tgg, logg, [Fe/H] and E(b-v) for 
every star within Kepler's field of view through a combi¬ 
nation of calibrated fluxes using {<7, r, f, filters similar 
to the original SDSS filters ([Fukugita et al. iri996fl and a 
narrow-band £>51 filter modeled after the Dunlap Obser¬ 
vatory £>£>51 filter. The catalog was originally used to 
inform target selection for the mission, but in the absence 
of a comprehensive spectroscopic survey of all ~ 150, 000 
Kepler stars, the catalog’s stellar parameters have been 
used in analyses of planet candidates. There are ongoing 
efforts to provide improved stellar parameters of Kepler 
targets by aggregating photometr y, spectroscopy, ast er- 
oseismology, and transit analyses (jHuber et al. 11201^ . 

The majority of false positive KPCs are expected to be 
caused by astrophysical sources rather than random or 
systematic errors, specihcally, EBs whose eclipse depths 
are similar to that ex pected from a transiting planet 
(jBorucki et al. Il2011br i. Fig. |T] demonstrates six of the 
most common sources of transiting KPC scenarios. In 
each panel, the larger (yellow) star is the suspected KPC 
host, and all objects within the panels are assumed to be 
within the aperture used to create the Kepler lightcurve. 
Each Kepler “optimal aperture” is variable, but is typ - 
ically many arcseconds in size ([Twicken et al. 1 120101) . 
The dashed circles represent a spectrograph fiber’s field- 
of-view (FoV, not to scale). The titles in each panel also 
denote, qualitatively, how often the given scenarios can 
be characterized by time-series RVs at modest precision 
(^ 100 m s“^ level). Note that in addition to stellar 
eclipses being diluted to look like giant planets, transits 
of giant planets can also be diluted to look like smaller 
planets. 

1 .2. Sourees Of False Positives 

Scenarios 1,2,4 and 6 all involve a physical companion 
orbiting the KPC host star. In these scenarios, RV obser¬ 
vations can detect the presence of grazing EBs (Scenario 
1), EBs that are diluted by light from a third star within 
the Kepler aperture, but resolved on-sky with the spec¬ 
trograph (Scenario 2), or consist of a very low-mass star 
(VLMS) or brown dwarf companion (Scenario 4), a ma¬ 
jority of the time. An important sub-category of Scenar¬ 
ios 1 and 4 include EBs whose orbits pr oduce only a sec¬ 
onda ry eclipse and no primary eclipse (iSanterne et al. I 
[201311 . These false positives may be more common for 
longer period KPCs, where a companion star in an ec¬ 
centric orbit is more likely to undergo secondary eclipse 
near periastron, but exhibit no primary eclipse. In addi¬ 
tion, the most massive, bona fide planets at short orbital 
periods will induce a Doppler velocity shift detectable at 
the ^ 100 m s“^ level (Scenario 6). In those rare cases, 
their planetary nature will be confirmed through the 
APOGEE RVs by phasing them to the Kepler-derived 
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Figure 1. The most common scenarios that can produce a 
lightcurve consistent with a transiting planet. The titles qualita¬ 
tively identify those scenarios that can be detected by a RV survey 
at the level of ~ 100 m s“^ (“Most”, “Some”, etc.). In each panel, 
the larger (yellow) star is the assumed KPC host star, while the 
dashed circles represent a spectrograph’s FoV (not to scale). All 
sources in each panel lie within the Kepler photometric aperture, 
which is typically many arcseconds in size. In Scenario 3, the term 
“unresolved” refers to the fact that the EB is unresolved in the 
Kepler aperture. In Scenario 6, only short-period, massive planets 
would be detected with APOGEE. Note that giant planets can also 
be diluted to look like smaller-sized planets in these scenarios. 


orbital period. 

Stellar systems that are either physical multiples or vi¬ 
sual companions with small separations on-sky, and con¬ 
sist of an EB, are represented by Scenario 5. In this 
scenario, the diluted EB is only detectable if the flux of 
at least one component of the EB pair is sufficiently high 
that it appears in the cross-correlation function, or if the 
combined mass of the EB pair induces a sufficient velocity 
shift on the third (brightest) star in the system. When 
the EB pair is composed of cooler K or M dwarf stars in 
the presence of a hotter primary, they are easier to detect 
in the NIR than in the optical, since the flux contrast is 
reduc ed in the H band (e.g., Kepler-16, iBender et al. I 
l2012t) . For a binary system composed of dwarf stars at a 
signal-to-noise of 100, secondaries with mass ratios down 
to ^ 0.1 are detectable in the H band (|Bender fc SimoiTI 
I2008D . while mass ratios are limited to ^ 0.5 in the op¬ 
tical. The detection limit for a given system depends 
on the number of stellar components within the aperture 
(e.g., is it a binary versus a triple system?), and whether 
any of those components are evolved (observing in the 
NIR is beneficial for components of differing ratios, 
not for brightness differences due to differing radii). 

Scenario 3 represents a ATepZer-unresolved EB, where 
the variable star is within the Kepler aperture, but is 
exterior to the spectrograph’s FoV relative to the KPC 
host star. This is the only scenario where RV observa¬ 
tions will be not be able to detect any false positives, 
unless the RV survey targets every star within a given 
KPC’s Kepler aperture. Fortunately, Scenarios 2, 3 and 
5 can sometimes be tested photometrically with time- 
series photom etry from the groun d at greater spatial res¬ 
olutions (e.g., IColon et al. Il201^ . In addition, these are 
also the scenarios that are more likely to be solved us¬ 
ing Kepler data alone, e.g., by searching for flux centroid 
shifts. 


1.3. Paper Outline 


In this paper, we introduce our program to ob¬ 
serve hundreds of KPCs using the Apache Point Ob- 
servatory Galact i c Evo lution Experime nt (APOGEE, 
Maiewski et al. I 120101) spectrograph (IWilson et al. 
gORl 12012 ! ) on the Sloan 2.5 m telescone (jGunn et al. 
20061) . recently finished as part of the Sloan D igital Sky 
Survey III (SDSS-III, lEisenstein et al. 1120111) and con¬ 
tinuing in SDSS-IV (2014-2020). Our program provides 
an efficient means of determining the false positive rate 
of KPCs due to physically-associated binary stellar sys¬ 
tems. At the same time, these spectra are used in a 
variety of projects concerning the false positives them¬ 
selves, including characterization of the orbits and mea¬ 
surements of the mass ratios for many of the spectro¬ 
scopic binaries (SBs), or orbital characterization of in¬ 
trinsically rare, massive {M > 10 Mjup), substellar com- 
pan ions such as brown dwarfs and massive gas giant plan¬ 
ets (jMarcv fc Butler II2000I l^hlmann et al. 1120 111) . For 
KPCs that remain viable, host star properties such as 
Teff, logg, and chemical abundances for dozens of ele¬ 
ments can be derived using the APOGEE spectra. 

In Section[2]we describe the APOGEE instrument and 
main survey, the methods used to derive RVs from its 
spectra, and its current RV precision floor. In Section 
ED we present RVs of five current and former KOIs ob¬ 
served during SDSS-III. Three of these happened to be 
observed as part of a separate APOGEE EB program, 
and we present some conclusions on the nature of those 
KOIs as a precursor to our larger KPC campaign. We 
also present the first results from our APOGEE-Kepler 
KOI campaign, using the (since confirmed) exoplanet 
host KIC 6448890 to test our long-term RV precision, 
and definitively identifying KIC 6867766 as a false posi¬ 
tive exoplanet. 

In Section 13.21 we compare the efficiency of a sur¬ 
vey using a high resolution, NIR, multi-object spec¬ 
trograph against other planet-hunting spectrographs: 
HARPS-North, which is a clone of HARPS-South with 
some improvement s (|Mavor et al. II2003D. Keck HIRE S 
(iVogt et al. I I1994D. SOPHI E ( Perruchot et al. I I2008D . 
and HET HRS ( Tull lITo^ . We demonstrate that by 
using a multiplexing instrument in the NIR to conduct 
a survey at modest RV precision (100 m s“^), false pos¬ 
itives can be identihed more efficiently compared to the 
single-object instruments, reserving telescope time on 
those other resources for confirmation of the remaining 
KPCs at significantly higher precision. In Section |4] we 
review other techniques for determining the false positive 
rate of Kepler KPCs, and highlight the science enabled 
by extracting abundances from the coadded APOGEE 
spectra. We summarize our findings in Section [5] 

2. APOGEE SURVEY OVERVIEW 

APOGEE is a survey of Milky Way stars using a multi¬ 
object, fiber-fed, NIR spectrograph housed in a vacuum 
cryostat, that can observe up to 300 objects simultane¬ 
ously, producing R ^ 22500 spectra covering a wave¬ 
length range of 1.51 — 1.68 fj,m using a volume phase holo¬ 
graphic grating mosaic. Detail s of th e inst rument design 
can be found in IWilson et al. I ()2010[ ) and IWilson et al. 1 
(l2012t) . Typically the instrument achieves a signal-to- 
noise ratio per pixel of 100 (AA ^ 0.1 — 0.17A) on an 
H = 11 star in a single visit (one hour of total in¬ 
tegration). Most stars are observed on a minimum of 

























































4 


Scott W. Fleming et al. 


APOGEE-KEPLER Fields In SDSS-III 



Figure 2. Kepler module footprints with APOGEE FoV (circles, 
2.98° diameters) from our SDSS-III Kepler EB and KOI programs. 
Five additional Kepler modules will be observed by APOGEE in 
SDSS-IV. 

three different nights, so that short-period binaries can 
be flagged. Each field on the sky is normally observed 
in multiples of three, ranging from 3 to 24 epochs, with 
brighter targets swapped for new stars after three obser¬ 
vations. Aluminum plug plates hold optical fibers that 
carry the star light from the telescope into the instru¬ 
ment. The primary science goal of the survey is to study 
the Milky Way by measuring radial velocities and chem¬ 
ical abundances of 10 ® red giant stars, but a vari¬ 
ety of additional science projects are included. A sum¬ 
mary o f the project can be found in lAllende Prieto et al. I 
(|2008bD . A detailed des c ription of the survey will ap¬ 
pear in iMaiewski et al. I (|2015D . Details of the target 
selection for th e surv ey in SDSS-III can be found in 
IZasowski et al. I (|2013D . 

The telescope’s field-of-view (FoV) covers a circular 
area 1.49° in radius, which matches well to the size of 
a given Kepler module. Fig. [5] shows the Kepler mod¬ 
ules’ footprints, along with the three SDSS FoVs for 
our programs observed during SDSS-IIL A total of 163 
KPC host stars are observed in the SDSS-III KOI field 
(blue). Those targets were selected from all KPCs that 
had H < 14, 153 of which were dispositioned as planet 
“candidates” as of August 2013 (four others were con¬ 
firmed exoplanets, six were not dispositioned yet). As 
can be seen, a single SDSS footprint covers most of a 
Kepler module’s FoV. In addition to the KPC hosts ob¬ 
served during SDSS-III, five additional Kepler modules 
will be observed in SDSS-IV. 

The APOGEE data p rocessing pipeline is described 
in iNidever et al. I (|2015l) . Basic steps include collaps¬ 
ing the detector exposures for each of the three NIR 
arrays from 3D data cubes to 2D images, flat fielding, 
aperture extraction, wavelength calibration, sky subtrac¬ 
tion, telluric correction, and measurement of RVs. The 
RVs currently calculated by the automated pipeline make 
use of a grid of synthetic spectr a calculated from AT- 
LAS9 stellar model atmospheres (jMeszaros et al. |l2012l : 
IZamora et al. 1120151) . The mean, internal RV precision 
of the pipeline-produced relative RVs is ^ 100 m s“^, al¬ 
though it does depend on signal-to-noise, spectral type, 
and level of residual systematics from the data process¬ 


ing. Critically, these pipeline-derived RVs only work for 
simple cross-correlation functions, and are expected to 
fail when there is contamination from multiple stellar 
spectra, as in the case for most binary stars. As such, we 
derive our own RVs using additional, interactive process¬ 
ing of the data. These steps include manually correcting 
residual OH sky emission lines, selecting templates from 
a finer grid, and interactively fitting cross-correlation 
peaks, which may often be asymmetric or have multi¬ 
ple components in the case of binary stars. We cal¬ 
culate uncertainties for our RV measurements f ollowing 
the m aximum-likelihood procedure laid out by iZucker I 
(|2003D . This approach derives an analytical relationship 
between the cross-correlation function and it’s first and 
second derivatives to account for uncertainty contribu¬ 
tions related to the sampling and sharpness of the corre¬ 
lation peak, and the signal-to-noise of the target and tem¬ 
plate spectra. The RVs are then fit u sing a custom wrap¬ 
per t o the RVLIN software package ([Wright &: Howardi 
1200911 ■ which includes the ability to fit both components 
of a double-line spectroscopic binary through an itera¬ 
tive approach, and forces some orbital parameters to be 
identical between both components (e.g., orbital period, 
eccentricity, epoch of periastron). In some cases we make 
use of our IDL-based Levenb erg-Marquardt fitting code 
used in iBender et al. I ([201211 . 

3. RESULTS 

3.1. Initial Case Studies 

We have selected five KOIs with diverse histories 
and current statuses to test and develop our analysis 
pipelines. Three of these targets were observed as part 
of an SDSS-III ancillary program studying Kepler EBs 
([Mahadevan et aL~l 120151 ) (hereafter MAH2015). These 
targets were at one point Kepler planet candidates, 
but were determined to be likely EBs by the time the 
MAH2015 observations began. Note that since these 
KOI hosts were observed through a different program, 
the total number of epochs for these targets is less than 
the number of epochs that the SDSS-Kepler KOI pro¬ 
gram obtains (the SDSS-III EB program obtained 3 — 6 
epochs for each target, compared to > 18 epochs for 
the SDSS-III KOI program). The two other KPC hosts 
presented here come from our SDSS-III Kepler KOI pro¬ 
gram. We summarize our findings on these targets to 
demonstrate the diversity of astrophysical configurations 
encountered in our spectroscopic observations. 

3.1.1. SDSS-III EB Program: KIC 1571511 

KIC 01571511 (KOI 362, Kp = 13.42, H = 

12.04), consists of an F -type dwarf and low-mass M 
dwarf ([Ofir et al. 1 1201211 , and was observed as part of 
MAH2015. The orbital period is 14.0224519 days and 
the radius as estimated in the NExScI KOI catalog is 
14.7 ± 6.4 i? 0 . The eclipses of such a low-mass star 
(~ 2 % decrease in flux during primary eclipse) are com¬ 
parable to those expected for a gas giant planet. Indeed, 
this star was origina lly suspected to be an overlooked 
gas g iant exoplanet (iCoughlin et al. I I201H lOfir et al~l 
l2012t) . In the specific case of KIC 01571511, there is a 
small secondary eclipse (^ 0.05%) detected in the Kepler 

^ http://exoplanetarchive.ipac.caltech.edu 
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Figure 3. APOGEE RVs of the Ae^ er EB KIC 01571511 and the 
best-fit model from lOfir et al.1 II2012I 1 shown as the solid line. A 
constant offset between the model and the APOGEE RVs has been 
applied to account for a zero-point offset. For the given Kepler pe¬ 
riod and epoch of transit, it is clear even with just three APOGEE 
RVs that the object is not a planet, because the change in RV over 
a short fraction of the orbit is much greater (~ 10 km s“^) than 
expected for a planetary mass. 


lightcurve, which can be used to derive an estimate of the 
relative Teff ratio between the primary and secondary, 
and can therefore be used to help determine whether the 
object is a likely stellar companion. However, there is 
no guarantee that an EB system with a primary eclipse 
will also show a secondary eclipse, nor that the secondary 
eclipse is detectable ev en with Kepler’s precision. In fact, 
iSanterne et al. I (j2012t l found that some of their false pos¬ 
itive KPCs were EBs in eccentric orbits for which only 
the shallower, secondary eclipses are present, but were 
mistaken as planetary transits across the primary. 

Fortunately, these EBs are fairly trivial to detect spec¬ 
troscopically, as demonstrated in Fig. |3l In this fig- 
ure , we plot t h e bes t-fit RV model from the analysis 
by lOfir et alTI (|2012h . noting that there is a typo in 
the value of loi in their Table 3 that is missing a mi¬ 
nus sign. We also plot the three APOGEE RVs obtained 
through the ancillary program (Table [1} . Only a con¬ 
stant offset between the model and APOGEE data is 
included to account for instrumental zero-point differ¬ 
ences. Even with three data points, the RV variation 
observed in the APOGEE RVs is inconsistent with a gi¬ 
ant planet, given the period and epoch of transit from 
the Kepler lightcurve, because the change in RV over a 
short fraction of the orbit is much greater (~ 10 km s“^) 
than expected for a planetary mass. These data also 
demonstrate that APOGEE is capable of producing RVs 
at the ~ 100 — 200 m s“^ for iJ ~ 12 stars based on the 
rms residual to th e well-determined orbital solution from 
10 fir et ain (|2012h . KIG 1571511 corresponds to Scenario 
4 in Fig. (TJ where a low-mass star generates an eclipse 
depth comparable to that expected from a giant planet. 

3.1.2. SDSS-III EB Program: KIC 3848972 

KIG 3848972 (KOI 1187, Kp = 14.49, H = 12.80) is 
listed in both the EB and KOI catalogs, and was observed 
as part of MAH2015. As a KOI, the target was listed as 
a “False Positive” in the Q1-Q8 catalog, but is currently 
absent in the Q1-Q16 catalog. The KOI Q1-Q8 catalog 


KOI Catalog, P = 0.370529, T, = 2454966.798190 



EB Catalog, P = 0.741057, = 2455001.998566 


-19.8 

-20.0 


Figure 4. Orbital phase-folded RVs for KIC 3848972, using the 
period and ephemeris values from the KOI Q1-Q8 Catalog (left) 
and EB Catalog (right). The period ambiguity arises from whether 
the system is treated as an object that only produces a primary 
transit (KOI solution), or an eclipsing binary that produces both a 
primary and secondary eclipse with similar depths (EB solution). 
No significant RV variation is seen beyond a few hundred m s“^, 
disfavoring a physically bound stellar companion at either of these 
orbital periods as the source of the Kepler signal. 

lists a period of P = 0.37052915 days, while the EB cat¬ 
alog lists a period that is twice as long (P = 0.741057 
days). The estimated radius reported in the KOI cata¬ 
log is 3.53 ± 0.93 P®. Where the EB Catalog assumes 
two nearly-equal eclipses from primary and secondary 
eclipse events, the KOI catalog reports half the orbital 
period and defines the secondary eclipse as undetected 
or absent. Multi-color, g r ound- based photometry was 
observed by iColon et al. I (I2012D using a tunable filter 
on the OSIRIS instrument on the 10.4-m Gran Telesco- 
pio Canarias (GTC). They find a consistent star-planet 
radius ratio (2a) in both their blue and red filters, but 
measure a statistically significant (5.8cr) difference in the 
eclipse depths. Interestingly, the color differences during 
eclipse suggest that the secondary component is bluer 
than the primary. 

Only three APOGEE spectra were obtained for this 
target as part of MAH2015, however, a check on bina- 
rity can still be performed provided the orbital phase 
coverage is reasonable. We conduct a one-dimensional 
cross-correlation using a K-type dwarf template (Peff = 
5000K). We do not see evidence for any significant ro¬ 
tational broadening greater than ~ 10 — 20 km s“^. Al¬ 
though the APOGEE spectra for this star are somewhat 
noisy, we find a single, very stable GCF peak with no 
RV variation greater than a few hundred ms“^ (Table 
[T|). There is no obvious correlation with orbital phase 
after folding on both the KOI and EB Gatalog periods 
and ephemerides, despite spanning ^ 80% of the KOI or¬ 
bital phase and ~ 20% of the EB Catalog orbital phase, 
respectively (Fig. |4]). If the signal was caused by a hot¬ 
ter (bluer) secondary orbiting a brighter primary, the 
expected RV amplitude should be many tens of km s“^. 

Another possible explanation is that the observed Ke¬ 
pler signal comes from an object transiting a low-mass, 
cool (red) star that is within the p hotometric aperture 
(constrained to be < 2" given the IColon et al. I (|2012f) 
GTC aperture), but too faint to detect spectroscopically 
with APOGEE in the presence of the bright primary 
star. In this case, the color-dependent transit depths 
are caused because the fainter, redder component is the 
one being transited, hence the overall color of the com¬ 
bined light appears to shift towards the blue during the 
transit event. Intrigued by this possibility, we obtained 
Keck adaptive optics (AO) imaging to search for a fainter 
companion that might be the source of the Kepler signal. 

The AO image was acquired on UT 2014 Jul 17, 
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Figure 5. Upper limits on secondary companion m ass from Keck 
AO im aging. The photometric aperture used by IColon et al. I 
()2Q12) is marked by the vertical line, and serves as the outer limit 
of where the eclipsing object might lie. 

using NIRC2 (instrument PI: Keith Matthews) and 
the Keck II Natural G uide Star (NGS) AO system 
(|Wizinowich et al.~l 1200011 . We used the narrow cam- 
era setting with a plate scale of 10 mas pixeP^, which 
provides a fine spatial sampling of the instrument point 
spread function (PSF). The observing conditions were 
excellent, with a seeing of 0".3. KIC 3848972 was ob¬ 
served at an airmass of 1.12. We used the Ks filter to 
acquire the image using a 3-point dither method. At each 
dither position, we took a total of 10 coadds composed 
of 5-second exposures. The total on-source integration 
time is therefore 150 seconds. 

The raw NIRC2 data were processed using standard 
techniques to replace bad pixels, flat-field, subtract ther¬ 
mal background, align, and coadd frames. We calculated 
the 5cr detection limit as follows. We first defined a series 
of concentric annuli centered on the star. For the con¬ 
centric annuli, we calculated the median and standard 
deviation of flux for pixels within these annuli. We define 
the 5cr detection limit as five times the standard devia¬ 
tion above the median flux. Representative 5tT detection 
limits are {1.6, 3.3, 5.0, 5.4} magnitudes for projected 
separations of {0".l, 0".2, 0".5, 1"}, respectively. 

We translate the 5a upper limits on companion bright¬ 
ness into upper lim its on companion mass using the SED 
models compiled in iKraus fc Hillenbrand I (|2007f ) , assum¬ 
ing the secondary is a bound companion and that differ¬ 
ential extinction between the two spectral types is mini¬ 
mal in the Kg band. For a given absolute Kg magnitude 
of the primary, the contrast curve from the Keck AO data 
gives a lower limit for the secondary’s absolute Kg mag- 
nitude, which we then interp olate into a mass using the 
iKraus &: Hillenbrand I (|2007l) models. We adopt primary 
spectral types of GO and K5 as conservative upper and 
lower limits based on our spectroscopic cross-correlation 
analysis. We are able to rule out any bound companions 
more massive than 0.2 Mq at the 5a level exterior to 0.2 
arcseconds (Fig. [5]). 

Alternatively, the transit signal could be caused 
by a fainter background or foreground star that is 
physically unassociated with KIC 3848972, but still 
within the Kepler photometric aperture. Following 


iMorton &: JohnsonH l)2011ll . we can estimate the prob¬ 
ability of having a blend source within a given aper¬ 
ture usiM a model of the galac tic population from TRI¬ 
LEG Al|3 dCirardOESZI tSi3)- We generate a TRI- 
LEGAL population in a one square degree area cen¬ 
tered on KIC 3848972 with the default settings (see 
Appendix We calculate a mean stellar density 
of 0.004581 stars arcsec“^ within Kepler magnitudes 
15.49 < Kp < 21.26. This magnitude range is cho¬ 
sen because it represents stars faint enough to be un¬ 
detected in the Kepler aperture but still able to pro¬ 
duce a transit depth of d = 0.00196. With this mean 
density, the probability of just finding a potential blend 
source within this magnitude range is 5.8% within 2"(the 
GTC photometric aperture) and just 0.36% within 0'.'5 
(the a rea least probed by t he Ky k AO images. Fig. [S]). 
From IMorton fc Johnson I (|2011ll . the probability that 
this blend source is an EB with a configuration that 
could mimic a transiting planet signal is on the order 
of 2.5E-4, so the probability of having a background EB 
as the source of this KPC is on the order of {9E-7, 1.4E- 
5} within (0.5, 2.0} arcseconds, respectively. Thus, the 
more likely EB false positive scenario is a bound eclipsing 
binary causing the transit signal. 

The lack of observed RV variability indicates this KOI 
is not due to a physically bound stellar companion orbit¬ 
ing the brightest component of the KIC 3848972 system, 
while the Keck AO images constrain any bound, diluted 
EB to be either within ^ 0.5" of the primary, or more 
than 5 magnitud es fainter than the p rimary in the Kg 
band. Given the IColon et al. I (j2012ri observations, we 
hypothesized that this KOI corresponded to Scenario 2 
or 5 in Fig. [U but our observations rule out Scenario 
2 and tightly constrain the separation of a diluted EB 
under Scenario 5. 

3.1.3. SDSS-III EB Program: KIC 3861595 

KIC 3861595 (KOI 4, Kp = 11.43, H = 10.27) is listed 
in both the EB and KOI catalogs, and was observed as 
part of MAH2015. As a KOI, the target was initially 
listed as a “Ealse Positive” in the Q1-Q8 catalog, and 
is currently listed as “Not Dispositioned” in the Q1-Q16 
catalog. The orbital period is 3.8493724 days and the 
estimated planet radius is 11.8 ± 1.6 i?®. Some ground- 
based observations have been conducted and reported 
at the Kepler Community Eollow-up Program (CFOP) 
websit({3- These inclu de several optical spectra fro m the 
TRES spectrograph (|Szentgvorgvi fc Fureszl [200^ that 
indicated the star was a rapid rotator (40-50 kms“^), 
and potentially variable at a level of a few hundreds of 
ms“^. Imaging from the 1-m Nickel telescope at Lick 
Observatory and Keck HIRES guider images show two 
nearby stars within ten arcseconds of the target. Both 
nearby stars appear to be approximately 6 magnitudes 
fainter than the target. 

In addition to six APOGEE spectra, MAH2015 ob¬ 
tained five optical spectra for this target using the High 
Resolution Spectrograph (HRS) on the Hobby-Eberly 
Telescope. The HRS was used in the 30,000 resolution 
mode, with the 316 g/mm grating at a central wave- 

^ http://stev.oapd.inaf.it/cgi-bin/trilegal_l.6 

^ https://cfop.ipac.caltech.edu/home/ 
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length of Ao = 5936 A. The HET spectra were re¬ 
duced using our optimal extraction pipeline described 
in MAH2015. We find a good template match (for both 
HET and APOGEE) using a mid-F spectral template 
rotationally broadened to 40 km s“^, in agreement with 
the CFOP notes from the TRES observations. The es¬ 
timated spectroscopic rotation rate of ^ 40 kms“^, 

comb ined with an estimated rotation rate of 5.65 -5.8 
days l)Hirano et al. |[20T^ [Rhodes fc Budding lIMT^ . re¬ 
sults in an equatorial radius of ^ 4.5 i?©, some¬ 

what larger than the spect roscopically determine d ra¬ 
dius of 2.727 ± 0.504 Rq (|Buchhave et al. II2012I1 and 
the Stellar Parameter Catalog’s value of 2.9921 q 743 Rq 
(jHuber et al. 1120141), which also uses the s pectroscopic 
stellar parameters of iBuchhave et ai~1 1)20121 1. but tie the 
stellar parameters to Dartmouth stellar evolution models 
()Dotter et al. l[2008ll . 

We cross-correlate the HET and APOGEE spectra 
with mid-F spectral templates rotationally broadened 
to 40 kms“^. For this target, we used subsections of 
the APOGEE spectrum (1.515-1.560 /rm, 1.586-1.605 
/rm, 1.615-1.635 /rm, and 1.6475-1.6775 /rm) to avoid 
several of the broadest lines. We find a single-peaked, 
broad cross-correlation function, with RV variation at 
the ^ 1 km s“^ level (Table |T]). We note that the 

RV scatter i s larger than m o st of the A and F stars 
observed by ILaerange et al. I (1200911 . However, phase¬ 
folding and fitting both sets of RVs at the orbital period 
and ephemerides found in the KOI and EB catalogs fails 
to resolve a signal of orbital motion from a bound com¬ 
panion at that period. 

In fact, we find that if we fit each set of RVs separately, 
fix the period and ephemeris to the KOI values, and force 
the eccentricity to zero, the best-fit RV semiamplitudes 
differ by a factor of 2.5 (Fig. IHl). A color-dependent 
semiamplitude may signal a blended spectrum (Scenar¬ 
ios 2 or 5); the redder component can affect the line 
shapes more significantly in the NIR, but such scenarios 
are particularly challenging to identify in rapidly rotating 
stars using only a handful of observations. Nevertheless, 
we undertake a full line bisector analysis of both HET 
and APOGEE spectra. After creating custom numerical 
stellar template masks for both the HET and APOGEE 
wavelength ranges, we calculate the bisectors of the cross- 
co rrelation funct i on, sim ilar to the procedure described 
in iWright et al. I (|20I3f) . We are limited to using three 
of the six APOGEE observations because they were the 
only ones observed on the same plug fiber, and there¬ 
fore should have the same intrinsic profile. The bisectors 
appear to be varying both in shape and position (indi¬ 
cating a cause other than bulk motion of the primary) 
and the bisector inverse slope (BIS) seems well-correlated 
with both RV and GCF FWHM (Fig. [7]). However, the 
rapidly-rotating nature of this star causes difficulty in 
establishing the GGF continuum, and complicates bisec¬ 
tor analysis. Any attempt to calculate errors on the BIS 
leads to overestimation, and therefore we are hesitant to 
quantify this result beyond saying that a blend scenario 
is possible. 

We can not definitively show the transit signal is 
caused by a spectroscopic blend with so few bisector mea¬ 
surements. Our HET RVs are consistent with a planetary 
companion in a circular orbit, but there are not enough 


KIC 3861595 - Comparison of APOGEE and HET RVs 



Phose (Wropped to Tq) 


Figure 6. Phase-folded RVs for APOGEE (red) and HET (blue). 
Each set were fit with the orbital period and transit ephemeris 
fixed to the KOI value. Eccentricity is forced to zero. We find 
that while both sets of RVs appear to be in-phase with the or¬ 
bital parameters, the RV semiamplitudes are quite different. This 
suggests the spectrum might be blended (Scenarios 2 or 5 in Fig. 
[T|l: the APOGEE spectra can be more sensitive to such a blend 
if the temperatures of the blended components are different. Line 
bisector variations also suggest a blend, but is not definitive due 
to the small number of observations. It is also possible that the 
uncertainties are underestimated due to the rapid rotation of the 
star (40 km s“^). 





Velocity (km/s) 



Velocity (km/s) 


Figure 7. Bisector analysis for the three APOGEE spectra that 
were on a common fiber (top row) and HET spectra (bottom row). 
The grey data point in the HET plots represents a low signal- 
to-noise observation. (Left: ) Bisectors of the cross-correlation 
function, shifted by measured radial velocities. Colors are based 
on bisector inverse slope (BIS) values. Percentage depth is a proxy 
for fiux, and does not span the full range from 0 to 1 because 
of continuum ambiguities. (Middle: ) Correlation between BIS 
and measured RV. (Right: ) Correlation between FWHM of the 
CCF, and measured RV. Uncertainties are determined from the 
measurement variation between echelle orders (for HET spectra) 
and between the three detectors (for APOGEE spectra). Note 
that difficulty in defining the CCF continuum probably leads to an 
overestimation of the BIS errors. 

RVs to make a firm claim. The APOGEE RVs contradict 
this claim, but RV uncertainties for rapid rotators have 
not been thoroughly vetted, so the APOGEE RV un¬ 
certainties reported in Table [I] coul d be underestimated. 
We al so note that an analysis by [Rhodes fc BuddinFI 
()20I4D found that an Algol-type background binary, ap¬ 
proximately 6.5 magnitudes fainter but within the Kepler 
photometric aperture, could produce a lightcurve similar 
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KIC 6448890, H = 10.33, RMS = 78.38 m/s 



Figure 8. The long-term RV rms of the confirmed exoplanet host 
star KIC 68448890 demonstrates we are able to achieve relative RV 
precision at the 100 m s“^ level. 

to what’s observed. This scenario, which corresponds to 
Scenario 3 in Fig. [l] also remains a possibility given the 
companions seen in the Lick and Keck images at this 
approximate flux ratio. This KOI is a prime example 
as to why it is sometimes necessary to obtain multiple 
spectra when searching for exoplanet false positives; ob¬ 
taining just a few spectra, even at orbital quadratures, 
may not be sufficient to confidently identify a blended 
stellar binary, especially for systems that are rapid rota¬ 
tors. Additional spectroscopic observations will be able 
to study the line bisectors and RVs in sufficient detail to 
determine the nature of this intriguing KOI. 

3.1.4. SDSS-III KOI Program: KIC 6448890 

KIC 6448890 (KOI 1241, Kp = 12.44, H = 10.33) is 
a system with two exoplanets that have been con firmed 
via transit timing variations (iSteffen et al. l[^013f) . The 
two planets have orbital periods of 10.5016 and 21.40239 
days, radii of 0.581 and 0.874 R.T„r.it,Br, a nd masses of 0.07 
and 0.57 Mj^piter (|Huber et al. 11201^ . The RV semi¬ 
amplitude is too small to be detectable with APOGEE, 
so this target {H = 10.33) is an opportunity to test the 
long-term RV precision level for our KOI program. Fig. 
[8] shows that the RV rms about the mean is 78 m s“^ 
over the entire baseline, in support of our stated goal to 
achieve a long-term (1-2 year), relative RV precision of 
~ 100 m s“^. The APOGEE RVs are reported in Table 

[II 

3.1.5. SDSS-III KOI Program: KIC 6867766 

KIC 6867766 (KOI 1798, Kp = 14.38, H = 12.99) 
was listed in the Q1-Q6 KOI catalog as an exoplanet 
candidate but has since been listed as “Not Disposi- 
tioned” in later catalogs. This KOI is also listed in the 
EB catalog. A shallow, 0.3% transit signal is present 
with no obvious secondary feature at an orbital period 
of 12.964725 days. The estimated radius from the KOI 
catalog is 9.65 ± 1.5 R^. This target was observed a to¬ 
tal of 25 times with APOGEE as part of our Kepler KOI 
program within SDSS-III. We found the best ID CCF 
template match using a Tgff = 5500 K, solar-metallicity 
dwarf rotationally broadened to 14 km s“^. Upon visual 
inspection, some of the CCFs were observed to be asym¬ 
metric, and in some cases evidence of a double-peaked 



Figure 9. Phase-folded ID CCF for KIC 6867766 as observed by 
APOGEE for a subset of the epochs. The CCF is clearly double- 
peaked, and the two CCF components have maximum separation 
near phases 0 ~ 0.25 and <f> ~ 0.75: a clear indication that the 
KOI planet signal is the result of a binary companion. 

CCF were present. In these situations, the parameters 
of our best ID CCF template are generally not reflective 
of any component in the system: a multi-dimensional 
cross-correlation analysis is required. 

To test whether the Kepler transits might be due to 
the binary showing up in the CCF, we calculate the ID 
CCFs for each APOGEE observation, normalize each 
CCF to a peak value of unity, and then sort them based 
on the orbital phase corresponding to each observation, 
phase-folding on the time of transit and orbital period 
as reported in the KOI catalog. We show only a few of 
the CCFs sampling different orbital phases in Fig. [S] for 
clarity. We find that this technique is quite effective at 
finding time-variable CCF changes indicative of blended 
SB2s. In the case of KIC 6867766, the double-peaked 
nature of the CCF is readily apparent, and phase-folds 
nicely with the KOI orbital solution (maximum separa¬ 
tion between CCF components near phases 0 0.25 and 

cj) 0.75 after phase-folding on KOI To, blended compo¬ 
nents near (j) 0 and 4> ~ 0.5). 

Upon further analysis, we determined that this sys¬ 
tem includes three stellar mass objects: a late F-type 
dwarf (KIC6867766A) with an early M-dwarf companion 
(KIC6867766B) in a ^ 13 day orbit, and a late G-type 
dwarf (KIG6867766G) with no discernible velocity mo¬ 
tion. We analyze d the APOGEE spe ctra with the TRI- 
COR algorithm (jZucker et al. I [19^ . which uses a 3D 
cross-correlation to derive the RV of three blended spec¬ 
tral components and their relative flux ratios. Our anal¬ 
ysis use d spectral ternplate s constructed from BT Settl 
models ([Allard et al. Il20flll . convolved to the APOGEE 
spectral resolution and sampling. We optimized the tem¬ 
plate stellar parameters to best match the observed spec¬ 
tra, using the peak correlation power from our highest 
S/N APOGEE spectra to assess the template match. All 
three components used templates with solar metallicity 
and \ogg = 4.5. The effective temperatures and rota¬ 
tional velocities were: A: 6200 K, 10 km s“^; B: 3500 K, 
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Figure 10. Three-component RV solution for KIC 6867766 using 
APOGEE RVs. Components A and B are an F+M EB pair that 
produces the Kepler transit signal, while Component C is a (likely) 
bound G dwarf companion with a long orbital period that dilutes 
the A-l-B eclipse. 

3 kms“^; C: 5200 K, 3 kms“^. We fixed the relative 
H band flux ratios to be C/A = 0.35 and B/A = 0.05, 
which are consistent with both the optimal flux ratios de¬ 
rived by TRICOR and the physical flux ratios expected 
for these stars. 

Table [1] contains the RVs we derived for each of the 
stellar components, while Table [2] lists the best (least- 
squares derived) spectroscopic orbital parameters for the 
13-day binary. We rejected some or all of the component 
RVs measured in five of the epochs due to peak pulling 
between the three components (generally A and C, but 
occasionally all three), and these are indicated in Ta¬ 
ble [T] with ellipses. Most of our APOGEE spectra have 
median S/N of 25 — 30 per pixel, but four (2014-04- 
11, 2014-04-14, 2014-05-11, 2014-06-18) had S/N < 10. 
We recovered the M-dwarf signal in each of these low 
S/N spectra and so retained them in our final orbital 
solution, albeit with larger RV uncertainties that reflect 
the poorer data quality. Fig. [10] shows the phase-folded 
RV curves for the AB binary, and residuals from the de¬ 
rived orbital solution. We also plot the measured RVs 
for K1C6867766C, which are flat and very close to the 
systemic velocity of the AB binary. This suggests that 
K1G6867766G is a bound companion to the AB binary 
in a very long-period orbit. We did not detect significant 
change in the RV of K1C6867766G, so cannot estimate 
its orbital period. We conclude that the K1G6867766 
is a likely hierarchical triple system, and an exoplanet 
false positive corresponding to Scenario 2 in Figure |T] 
(although in this case, the tertiary star is also within the 
spectrograph’s FoV). 

3.2. Survey Efficiency 

Several high-precision RV instruments are available in 
the northern hemisphere to confirm Kepler exoplanet 
candidates by achieving RV precisions of a few ms“^, 
including HARPS-North, Keck HIRES, SOPHIE, and 
HET HRS. These instruments are clearly capable of con¬ 


ducting a reconnaissance survey of Kepler KPCs for false 
positives, but their telescope time is most effectively 
spent observing robust exoplanet candidates for con¬ 
firmation and characterization purposes. Smaller tele¬ 
scopes have been used to measure RVs for hundreds of 
KOIs through the Kepler GFOP (e.g., McDonald 2.7m, 
Tillinghast 1.5m), and while they have helped identify 
the best candidates for higher precision RV follow-up, 
they are limited to a single target at a time: most KOIs 
have just one or two observations from these facilities. To 
compare the efficiency of APOGEE against the larger- 
aperture telescopes mentioned above, we calculate the 
integration time per target required to achieve a photon- 
limited RV precision of 100 m s~^. A qua lity factor Q is 
calculated following iBouchv et al. I (| 2001 [ ). adapting in¬ 
strument parameters summarized in Table |3| The qual¬ 
ity factor represents the fundamental RV information 
content of a spectrum, which depends on the number, 
depths, and widths of spectral features. The Q values 
are calculated using BT Settl stellar models for a range 
of stellar effective temperatures (Teff), adopting a surface 
gravity logg = 5.0, a solar metallicity, and no rotational 
broadening. Including a rotational broadening will af¬ 
fect higher resolution instruments the most, resulting in 
lower Q values. 

From these Q values, it is then possible to calculate 
the required integration time per object (an “effective 
integration time”) to achieve a given, photon-limited RV 
precision (here taken to be 100 m s“^), via: 

^targets 

where t is the effective integration time, c is the speed 
of light, (Trv is the desired RV precision, Q is the qual¬ 
ity factor, R is the telescope’s effective aperture radius, 
e is the total throughput (as a percentage) of the tele¬ 
scope and instrument, F is the flux in photons per second 
per unit area, tover is the overhead per integration, and 
n-targets is the number of targets observed per integration. 
Since we are operating in the photon-limited case, we do 
not include readout noise for the instruments, nor con¬ 
sider sources of systematic uncertainties such as residual 
moonlight contamination (worse in the optical), or resid¬ 
ual telluric lines and sky emission lines (worse in the NIR) 
- our interest is in calculating the photon limited case for 
a direct comparison. 

The fluxes are calculated from the BT Settl models 
after convolving with the appropriate filter transmis¬ 
sion function: Johnson/Bessell 1^ for HARPS-North, 
Keck HIRES, SOPHIE, and HET HRS; 2MASS for 
APOGEE. The model fluxes are scaled to the zero mag¬ 
nitude level usi ng zero-level f l uxes f rom iBessell et al. I 
(|I998l) in V and ICohen et al. I (I2003D in 2MASS, which 
are then further scaled to a desired apparent magnitude. 
We calculate the effective integration time for each in¬ 
strument by including an estimated overhead for detector 
readout, telescope slew, target acquisition and calibra¬ 
tions. Since overhead times between integrations depend 

® http://www.ctio.noao.edu/~points/SOIFILTERS/fliters/ 
maintext.html 

® http://www.ipac.caltech.edu/2inass/releases/allsky/doc/ 
sec3_lbl.html 
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For H = 13.5 




Figure 11. Total time (integration + overhead) required to 
achieve a photon-limited RV precision of 100 m on a per target 
basis, as a function of stellar Tgff. Note that at this scale, HARPS- 
North and Keck HIRES are nearly indistinguishable. APOGEE’S 
multiplexing capability is the driving factor in reducing the time 
per target here. 

on a variety of factors, we adopt three minutes as an 
average overhead time for HARPS-North, Keck HIRES, 
SOPHIE and HET HRS. For APOGEE, we use a mini¬ 
mum integration time of 66.7 minutes and an overhead 
time of 15 minutes, equivalent to the current survey’s 
integration time per visit for each field and its average 
overhead time between fields. We do not reduce the in¬ 
tegration time for APOGEE below this minimum value 
because the total integration time is not based on any 
one target’s brightness, but rather an overall integration 
time required for all targets in a field. 

We calculate the required integration times t at a spe¬ 
cific H (and corresponding V) magnitude for both the 
single object and multi-object instruments, to allow for 
more direct comparisons. We use the median KPC host 
star flux level of H ~ 13.5, along with the appropriately 
scaled V fluxes for the optical instruments. The average 
number of KPG host stars within each APOGEE field 
that have H < 13.5 is 89, so we use this as the number 
of KPG hosts that can be observed simultaneously with 
apogee’s multiplexing capability for the purposes of 
comparing against the single-object instruments. All of 
these input parameters are summarized in Table [3l 

Eig. [TT] displays the required integration times per tar¬ 
get to achieve a photon-limited RV precision of 100 m s“^ 
for a variety of spectral types. At this RV precision, 
APOGEE is approximately three times as efficient on 
a per target basis, primarily due to APOGEE’S abil¬ 
ity to observe multiple KPG hosts simultaneously, and 
the fact that at this RV precision, the other instruments 
are dominated by overheads per target. As the targeted 
precision level increases, the other instruments become 
increasingly more efficient compared to APOGEE, re¬ 
flective of the fact that they are no longer dominated by 
overheads. These basic calculations serve to demonstrate 
how the multiplexing capability of APOGEE enables an 
efficient survey at modest RV precision compared to the 
single-object instruments, and that the telescope time 
for those instruments is best spent on achieving high RV 
precision to confirm new exoplanets. 

Eig. [T2]plots the orbital period versus H magnitude for 


Figure 12. Orbital period versus H magnitude for the current 
KPCs. The dashed, red lines demarcate H = lA and P = 100 days. 
Host stars brighter than this H limit and candidates with periods 
below 100 days represent > 80% of the current KPC catalog, and 
are particularly well-suited for APOGEE to characterize the orbits 
of binaries. Our SDSS-III and SDSS-IV KOI programs observe 
targets over more than one year baselines, however, and can detect 
stellar companions out to the longest KOI orbital periods. 

the current KPC catalog. Host stars with H < lA and 
P < 100 days are particularly well-suited for APOGEE 
to characterize any binary star orbits, and represent more 
than 80% of the current KPCs. Note, however, that even 
KPCs with longer orbital periods, extending to at least 
a few hundred days, can be identified as binaries even if 
the observing baseline does not cover the entire orbital 
period. If a more conservative limit of 77 < 13 is applied, 
more than 47% of the KPCs would be included. Scenario 
5 in Fig. [T] relies most heavily on achieving high signal- 
to-noise ratio observations, and thus has the brightest 
limiting magnitude within a survey; however, it is also 
the false positive scenario for which APOGEE is least 
sensitive. 

4. DISCUSSION AND COMPARISON TO PRIOR WORK 
4.1. Comparison With Other Efforts 

A variety of techniques exist to determine whether a 
given KPC might be a false positive using only Kepler 
photometric information. Examination of the transits 
can be done to ensure that the odd/even transits have 
the same depth, that there are no ellipsoidal variations, 
and that the positions of the flux centroids do not vary 
with brightness changes, all of which are signs that the 
trans i t event may be due to a diluted EB (|Steffen et al. I 
l201Ct [Torres et al. Il201ll) . Another method is to gen¬ 
erate a large grid of synthetic EB blend scenarios and 
compare the models to the ob served lightcurves v ia a 
analysis (e.g., BLENDER, iTorres et, al. Il2011^ . al¬ 
though the technique can be computationally expensive 
and difficult to apply to all KPCs en masse. In addition, 
imag ing surveys (e.g., iHowell et al. l[20TltlAdams et al. I 
12012 ) can be used to inform the photometric analyses de¬ 
scribed above, particularly to identify fainter stars that 
exist within the lightcurve aperture, and to search for 
wide stellar companions to study any relationships be- 
tween exoplanet p roperties and host star multiplicity 
(|Wang et al. I[20^ . 

Another technique makes use of stellar population 
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synthesis and Galactic models to estimate the prob- 
ability that a given tran sit signal is due to an EB. 
iMorton fc Johnsonl l)2011h used such a technique to esti¬ 
mate the false positive probability (FPP) of KPGs, and 
concluded that 90% of the KPGs had FPPs < 10%, a 
result that was used by other authors in subsequent sta¬ 
tistical analyses of the KP G planet candidates. Howe ver, 
the technique presented in iMorton fc: Johnson I (1201111 re¬ 
lied on KPGs to be vetted to the fullest extent possible 
using the Kepler photometry, specifically, the removal 
of V-shaped transits and searches for faint secondary 
eclipses, whereas the table of KPGs presented in their 
paper was not limited to these pre-vette d KPGs. 

An updated version of the technique bv lMorton~l 1)201211 
accounted for more false positive scenarios and clarified 
the importance of pre-vetting KPGs before performing 
the statistical analysis. The updated implementation 
can be run fairly quickly (of order 10 minutes per star), 
and has been verified by testing it on confirmed KPGs 
and known false positive KPGs; however, the technique 
works best with additional observations (imaging to de¬ 
tect close, visual companions and at least one high reso¬ 
lution spectrum to get coarse stellar parameters). In ad¬ 
dition, the Bayesian modeling is dependent on a variety 
of model assumptions regarding Galactic structure, stel¬ 
lar population synthesis, distribution of binary proper¬ 
ties, and the frequency of exoplanets for various types of 
stars. While the framework explicitly accounts for such 
assumptions through the adopted priors, and is fairly 
trivial to update when new knowledge is obtained about 
any of these distributions, direct spectroscopic or photo¬ 
metric observations of false positives are the least model- 
dependent approach to derive the false positive statistics 
of KPGs. 

The Kepler CFOP program has been conducting 
spectroscopic and imaging campaigns to identify blend 
sources of KOIs, and, as previously mentioned, have col¬ 
lected thousands of RV measurements for hundreds of 
KOIs, in addition to high-resolution imaging to search for 
faint companions unresolved in the Kepler photometric 
aperture. Our program selected targets independently 
of those observed by the GFOP program, since our goal 
was to minimize selection bias in the KOIs we observe. 
Our program complements the CFOP in a variety of 
ways. The CFOP imaging data can identify wide stellar 
companions, while our RV measurements over 18-month 
baselines can detect linear RV trends of intermediate- 
period companions that are unresolved by AO or lucky 
imaging. The abundances of the ~15 elements within 
apogee’s H band spectral window can be compared 
and contrasted with elements accessible in the CFOP 
optical spectra. 

Many of the CFOP spectroscopic observations have ob¬ 
tained just one or two spectra near predicted quadrature 
phases, where any contaminating spectral lines from a 
stellar companion are at maximum separation. This is of¬ 
ten sufficient to identify a subset of false positive scenar¬ 
ios: eclipsing stellar companions orbiting the KIC star, 
for example. In contrast, each of our KOI targets in 
SDSS-III/IV (except those few that were part of our Ke¬ 
pler EB program) obtain more than twenty APOGEE 
RVs, which allow us to fully characterize the orbits of 
any bound companions causing the false positive transit 
signals. In the case of eclipsing low-mass secondaries. 


this further enables a study of the fundamental mass- 
radius relationship for K and M dwarfs, since precise 
mass ratios at the 1% level are achievable. Our mul¬ 
tiple RV measurements can also be used to search for 
(and place limits on) the presence of any longer-period, 
non-transiting companions (stellar or otherwise), for the 
study of multiplicity amongst Kepler planet hosts. 

Studies have found that the false positive rates for 
various subsets of KPGs are larger than the ones found 
from the IMorton fc Johnson I (|2011h study, whose quoted 
statistics are only valid for fully pre-vetted KPGs, and 
that this rate might differ dependi ng on the orbi t al pe - 
riod and transit depth of the KPGs. IColon et al. I (120121 ) 
made use of multi-color differential photometry to test 
for false positives due to diluted EBs whose components 
have sufficiently different colors. They observed a total 
of four KPGs that had short periods (P < 6 days) and 
small radii {Rp < 5 P©), and found evidence that two 
of the four were likely due to diluted EBs, excluding an 
overall false positive r ate of 10% with 99% confidence. 
iSanterne et al. I (I2012D collected spectroscopic RVs of 33 
giant planet KPGs using the SOPHIE spectrograph and 
found a false positive probability of at least 35% within 
their sample, where a majority of false positives were due 
to EBs. Their sample size of 33 KPGs was partly limited 
by their telescope resource: a single-object spectrograph 
observing in the optical using a 2m-class telescope, cor¬ 
responding to an effective magnitude limit in the Kepler 
bandpass of Kp < 14.7. This magnitude limit removes 
almost half of the total KPGs. Utilizing a multi-object, 
NIR spectrograph, such as APOGEE, increases the rate 
of data collection while also increasing the total number 
of KPGs able to be observed. 

4.2. Abundances of KOI Host Stars 

Beyond identifying false positive KPGs as binaries, 
a variety of additional science projects can be done 
with the NIR APOGEE spectra. One such example 
is the study of chemical abundance patterns in planet 
host stars compared to stars not known to host plan¬ 
ets. One of apogee’s primary goals is to measure the 
chemical abundances of many elements to study stel¬ 
lar populations within the Milky Way. These abun¬ 
dances are measured using the APOGEE Stellar Pa- 
rameters and Chemical Ab undances Pipeline (ASPCAP, 
iGarcia Perez et al. I 1201511 . which consists of a suite 
of software codes to analyze, in an automated fash¬ 
ion, the APOGEE spectra. The main component of 
the code is FERREQ, a Fortran optimization code that 
searches for the set of parameters that best match 
each APOGEE spectrum. FERRE was originally devel- 
oped in the context of low-resolution SPSS spectro scopy 
(jAllende Prieto 1 12004 lAllende Prieto et al. 11200^ , and 
has s ubsequently evolved and been used in other con¬ 
texts (lAllende Prieto et al. 1 1200^ 120091: iBrown et al. I 
[Ml iKilic et al. 1 1201211 . The APOGEE band (1.5-1.7 
^m) is rich in transitions from many elements in cool 
stars. Abundances for 15 elements can be derived from 
sufficiently high resolution (R > 20000) and signal-to- 
noise per pixel (S/N> 100) spectra in this spectral win¬ 
dow: C, N, O, Na, Mg, Al, Si, K, Ca, Ti, V, Mn, Fe, Ni, 
Cr. A S/N level approaching 100 per pixel is expected 

^ FERRE is available at http://hebe.as.utexas.edu/ferre 
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ASPCAP [M/H] vs. Buchhave et al. (2014) 



Buchhove et ol. (2014) [M/H] 

Figure 13. Compariso n of [M/H] values from ASPCAP DR12 and 
IBuchhave et al.~1 1120141 ) for stars in common. The ASPCAP values 
are taken from the "’uncalibrated” (initial fit) values (“FPARAM” 
array). We find that a majorit y of these targets fall within 0.1 dex 
of the IBuchhaye et al. I II2014I ) values (grey dashed lines). These 
boundaries also represent the typical la uncertainties for both sets 
of values. There are not enough targets with “externally cali¬ 
brated” yalues (“PARAM” array) in ASPCAP DR12 to make a 
statement about how those compare to IBuchhave et al. I II2014I ). 

to be achieved out to i? ~ 14 by coadding the multiple 
visits. 

A detailed analysis of elemental abundances for KOls 
with and without exoplanets is beyond the scope of this 
introductory paper. To get some sense of the metal- 
licites coming from the automated ASPCAP pipeline, 
we have compared the “uncalibrated” [ M/H] values from 
ASPC/ AP in DR12 with those found in IBuchhave et alTI 
(j2014l) for stars in common. The “uncalibrated” values 
are the parameters that come from the initial fit. A 
subset of targets have additional, external calibrations 
applied to their stellar parameters (such as metallici- 
ties of clusters from the literature). We refer the reader 
to the ASPCAP DR12 documentatioifl for full details. 
There are not enough targets that have had external 
calibrations applied to make any statement regarding 
their agreement. However, we find a total of 128 KOI 
host stars observed in common that have “uncalibrated” 
[M/H] values in DR12. We find the agreement to be 
promising (Fig. (TS]): a majority of targets agree within 
0.1 dex (grey lines in Fig. [TS]), despite the fact that AS¬ 
PCAP has been calibrated primarily to work on bright 
giants. We have not made any cuts based on ASPCAP 
processing flags; this is a compa rison using DR12 [M/H ] 
values “as they are” versus the IBuchhave et al. I (|2014D 
values. As such, the relation in Fig. [13] should be consid¬ 
ered as preliminary, and likely to be improved upon in 
future analyses. 

There are tantalizing hints for different heavy ele¬ 
ment patterns in planet hosts relative to the field that 
could be induced by preferential removal of heavy ele¬ 
ments in the disk. It’s been shown that stars hosting 
Jovian-mass planets tend to be more me tal-rich than 
stars with only Neptunian-mass planets (iSousa et al. 1 
120081 : IChezzi et al. ll20I0HSousa et al. II20IID . The over- 
all shift in values of [Fe/H] between stars with Jovian- 


http://www.sdss.org/drl2/irspec/parameters/ 


versus Neptunian-mass planets is 0.20 dex, which is sig¬ 
nificant and indicates that the metallicity populations for 
stars with Jovian-mass planets are not the same as those 
which host the smaller Neptunian-mass planets. Oth¬ 
ers have used the Kepler sample to extend the planet- 
metallicity correlation down to terrestrial-sized planets 
([Buchhave et al. 1l2012t) , and have found that terrestrial¬ 
sized pla nets fall into well-defin ed host-star metallicity 
regimes ([Buchhave et al. II2014D . These results suggest 
that metallicity may also influence the distribution of 
planetary masses within extrasolar systems. In addi¬ 
tion to stellar metallicity, there are also suggestions that 
stellar mass plays a role, such that the dominant plane- 
tary mass decreases a s the parent star’s m ass decreases 
([Cirezzi et al. I I20I0I: iMavor et al. 1 120111) . at least for 
main sequence stars. While massive subgiant and giant 
stars show trends similar to main se quence stars, low- 
mass giants show a different behavior ([Maldonado et aL~l 
120131) . It is likely that some combination of stellar mass 
and metallicity influences the type of planetary system 
that will form (iJohnson et al. 1120101) . APOGEE can 
provide a statistically significant control sample for such 
studies. 

In addition to overall metallicity and stellar mass play¬ 
ing a role, the detailed chemistry of the parent cloud in 
which the system forms may also hold clues to further 
understanding planet formation. Recent findings sug¬ 
gest that specific abundance patterns, such as Mg/Fe, 
may influence the likelih ood that a star hosts an u nder- 
lying planetary system ([Adibekvan et al. 1l2012aD . and 
that enhancement in alpha elements may favor the for¬ 
mation of ro cky planets, even for st ars with low iron 
abundances (lAdibekvan et al. 1l2012bD . The C/0 ratio 
in the parent cloud is also found in some studies to be en- 


riched in planet hosting syst e ms, with C/0 ratios > 0. 8 
([Delgado Mena et al. 1 120101 : iPetigura fc Marcy~l 120111). 
howey er, these results haye been questioned by iFortney I 
"Ml). and h a ye not been confir me d by other groups 


( Nissen ] 120131 : iTeske et al. 1 12014D . iBrugamyer et al. 1 


201 If ) find that the silicon abundance (and not the 


oxygen abundance) is a key element, as they find that 
their planet detection rate depends strongly on the sil¬ 
icon abundance of the host star. A difference in the 
Si abundance is also found for the XO-2 binary host 
stars, whe re XO-2N is found to be enhanced r elatiy e 
to XO-2S (iTeske et al. Il2 015). [Melendez et al. 1 (l2 009f). 


iConzalez Hernandez et al. 1 ( 20101 ). and LSchuler et al. 1 

(120111) present intriguing results suggesting that low- 
amplitude chemical signatures point to selectiye accre¬ 
tion or depletion of refractory and volatile elements in 
stellar exoplanetary hosts. In particular, trends of abun¬ 
dances with condensation temperature (Tcond) are used 
as diagnostics, and these can be defined from the abun¬ 
dances of the 15 chemical elements covered by APOGEE, 
which include C, N, O (volatiles) and Si, Ti and Al (re¬ 
fractories) . The investigation of such trends in samples of 
Kepler stars with confirmed planets of different masses, 
and including the smallest planets to date, provides an 
unprecedented database in order to probe the importance 
of Tcond trends in this context. 


5. SUMMARY 

In this paper, we highlight the importance of an RV 
survey of KPCs to better determine the false positive 
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rate, and demonstrate that APOGEE can efficiently con¬ 
duct such a survey of KPCs to identify, and in many 
cases characterize the orbits of, false positive KPGs. We 
have shown that the APOGEE instrument is capable of 
achieving an RV precision of ^ 100 m s“^ using observa¬ 
tions of the confirmed exoplanet host star KIC 6448890, 
as well as the Kepler EB KIC 01571511, which produces 
planet-sized eclipses and has an H magnitude similar to 
many KPCs. We find that the transit signal of KIC 
3848972 is not caused by a blue, stellar secondary orbit¬ 
ing the primary star, and do not find any evidence of a 
faint, red companion in Keck AO images that could be 
the source of the Kepler transit events. Further inves¬ 
tigation is merited before the true nature of this KOI 
can be confidently identified. We find HET RV varia¬ 
tions that phase to the KOI period and ephemeris for 
KIC 3861595, but our APOGEE RVs are inconsistent 
with the HET RV semiamplitude, and we find evidence 
of line bisector variations. This target was part of our 
EB program in SDSS, so we only have a few spectra 
to work with, and the RV uncertainties have not been 
fully vetted for rapdily rotating stars such as this one. 
As such, we can not definitively determine the nature of 
this KPC, and urge more spectra be obtained to exam¬ 
ine both the RV and line bisector variations. Finally, we 
find that KIC 6867766 is a triple system, composed of 
an F-l-M EB and a wide, bound G dwarf tertiary. The 
F-l-M EB phases to the KOI period and ephemeris, and 
the diluted eclipses are the source of the KPC transits. 
As such, we can confidently identify this KOI as a false 
positive exoplanet candidate. 

Not only can the data from such a survey be used to de¬ 
termine the false positive rate of KPCs and vet the sam¬ 
ple to identify the best candidates for high-precision RV 
observations, but it will enable ancillary science projects 
in fundamental stellar astrophysics though observations 
of EBs, studies of intrinsically rare short-period com¬ 
panions (such as brown dwarfs), and detailed chemical 
abundances of exoplanet host stars. At the precision 
level of 100 m s“^, APOGEE is a more efficient instru¬ 
ment compared to HARPS-North, Keck HIRES, and SO¬ 
PHIE, due to its multiplexing capability and because the 
single-object spectrographs are dominated by overheads. 
Our survey to detect false positives refines the target 
selection for higher precision RV instruments, enabling 
them to focus on the best exoplanet candidates. It will 
allow for improved statistical studies of the Kepler exo¬ 
planet population by determining the false positive rate 
of KPCs due to physically-bound binaries, as well as any 
trends in the false positive rate with orbital period or 
stellar properties. 

We thank the anonymous referee for detailed com¬ 
ments that improved the quality of this publication, es¬ 
pecially their suggestion to conduct more analysis on 
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APPENDIX 

TRILEGAL INPUT PARAMETERS 

The parameters used to calculate the background EB blend probability of KIC 3861595 are summarized in Table SI 
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Table 1 

RVs for KIC Stars - All RVs in km s”! 


KIC ID 

BJD.TDB 

RVa 

la 

RVs 

la 

RVc 

la 

Instrument 

1571511 

2455811.61304 

-24.401 

0.153 





APOGEE 

1571511 

2455840.59327 

-26.348 

0.115 





APOGEE 

1571511 

2455851.57845 

-18.927 

0.105 





APOGEE 

3848972 

2455811.61297 

-19.943 

0.157 





APOGEE 

3848972 

2455840.59327 

-20.161 

0.153 





APOGEE 

3848972 

2455851.57848 

-19.641 

0.150 





APOGEE 

3861595 

2455789.84195 

-23.052 

0.091 





HET 

3861595 

2455796.82710 

-23.827 

0.075 





HET 

3861595 

2455797.80608 

-23.204 

0.081 





HET 

3861595 

2455801.80843 

-23.194 

0.098 





HET 

3861595 

2455803.80347 

-24.285 

0.103 





HET 

3861595 

2455813.70317 

-22.543 

0.566 





APOGEE 

3861595 

2455823.72718 

-23.628 

0.469 





APOGEE 

3861595 

2455840.66180 

-22.930 

0.541 





APOGEE 

3861595 

2455849.57900 

-24.600 

0.434 





APOGEE 

3861595 

2455851.64939 

-21.337 

0.433 





APOGEE 

3861595 

2455866.56998 

-21.608 

0.461 





APOGEE 

6448890 

2456368.99828 

-55.391 

0.117 





APOGEE 

6448890 

2456411.92027 

-55.513 

0.105 





APOGEE 

6448890 

2456557.73343 

-55.598 

0.103 





APOGEE 

6448890 

2456559.72336 

-55.644 

0.106 





APOGEE 

6448890 

2456560.72108 

-55.571 

0.104 





APOGEE 

6448890 

2456584.63225 

-55.582 

0.104 





APOGEE 

6448890 

2456585.63076 

-55.644 

0.105 





APOGEE 

6448890 

2456757.89294 

-55.622 

0.107 





APOGEE 

6448890 

2456758.90229 

-55.801 

0.142 





APOGEE 

6448890 

2456760.90571 

-55.586 

0.112 





APOGEE 

6448890 

2456761.87281 

-55.573 

0.139 





APOGEE 

6448890 

2456762.86860 

-55.621 

0.111 





APOGEE 

6448890 

2456763.88112 

-55.581 

0.109 





APOGEE 

6448890 

2456783.83567 

-55.712 

0.112 





APOGEE 

6448890 

2456784.82195 

-55.781 

0.133 





APOGEE 

6448890 

2456785.82543 

-55.702 

0.108 





APOGEE 

6448890 

2456786.79845 

-55.590 

0.113 





APOGEE 

6448890 

2456787.80934 

-55.640 

0.107 





APOGEE 

6448890 

2456788.84307 

-55.679 

0.118 





APOGEE 

6448890 

2456812.74509 

-55.620 

0.111 





APOGEE 

6448890 

2456814.75547 

-55.615 

0.114 





APOGEE 

6448890 

2456815.78552 

-55.607 

0.107 





APOGEE 

6448890 

2456816.76627 

-55.710 

0.119 





APOGEE 

6448890 

2456817.76198 

-55.632 

0.109 





APOGEE 

6448890 

2456818.76458 

-55.609 

0.110 





APOGEE 

6448890 

2456819.76222 

-55.567 

0.109 





APOGEE 

6448890 

2456820.75601 

-55.666 

0.108 





APOGEE 

6867766 

2456557.73337 







APOGEE 

6867766 

2456559.72331 

38.598 

0.383 

-60.127 

2.882 

6.928 

0.664 

APOGEE 

6867766 

2456560.72103 

40.497 

0.365 

-65.538 

6.014 

8.042 

0.653 

APOGEE 

6867766 

2456584.63222 

28.939 

0.399 

-36.948 

3.811 

7.788 

0.780 

APOGEE 

6867766 

2456585.63072 

37.993 

0.446 

-64.834 

4.138 

8.956 

0.754 

APOGEE 

6867766 

2456757.89298 







APOGEE 

6867766 

2456758.90233 

1.277 

1.679 

46.983 

9.919 

8.567 

3.754 

APOGEE 

6867766 

2456760.90575 

-20.442 

0.482 

88.944 

3.869 

10.516 

0.778 

APOGEE 

6867766 

2456761.87284 

-17.315 

1.389 

79.738 

11.967 

10.500 

3.354 

APOGEE 

6867766 

2456762.86864 

-10.374 

0.432 

66.784 

3.542 

9.447 

0.760 

APOGEE 

6867766 

2456763.88116 



28.563 

3.463 



APOGEE 

6867766 

2456783.83568 

15.664 

0.590 

-3.902 

3.310 

11.547 

1.018 

APOGEE 

6867766 

2456784.82197 



37.085 

3.697 



APOGEE 

6867766 

2456785.82544 

-12.960 

0.422 

73.438 

3.067 

9.843 

0.759 

APOGEE 

6867766 

2456786.79846 

-21.181 

0.542 

89.949 

5.130 

9.443 

0.887 

APOGEE 

6867766 

2456787.80935 

-20.018 

0.398 

92.288 

4.554 

9.495 

0.656 

APOGEE 

6867766 

2456788.84309 

-12.259 

1.175 

55.283 

9.284 

7.295 

1.769 

APOGEE 

6867766 

2456812.74507 

-20.349 

0.501 

90.837 

6.466 

8.471 

0.870 

APOGEE 

6867766 

2456814.75546 

-10.567 

0.455 

70.185 

6.978 

9.955 

0.891 

APOGEE 

6867766 

2456815.78551 



32.014 

3.443 



APOGEE 

6867766 

2456816.76626 

16.157 

0.886 

-7.116 

4.691 

9.461 

1.656 

APOGEE 

6867766 

2456817.76197 

27.340 

0.445 

-35.030 

4.168 

11.971 

1.059 

APOGEE 

6867766 

2456818.76456 

37.006 

0.398 

-59.225 

3.500 

8.936 

0.702 

APOGEE 

6867766 

2456819.76220 

40.742 

0.422 

-65.495 

3.783 

9.396 

0.780 

APOGEE 

6867766 

2456820.75599 

38.303 

0.356 

-57.968 

3.165 

8.577 

0.668 

APOGEE 
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Table 2 

KIC 6867766 A+B Orbital Parameters 


Parameter 

Value 

ICT 

P (days) 

12.964712 

(fixed at Kepler value) 

Tp 

2456746.58 

0.21 

e 

0.0553 

0.0054 

U} (deg) 

128.6 

5.7 

(km s~^) 

30.77 

0.14 

Kq (km s“^) 

79.4 

1.2 

7 (km s~^) 

10.87 

0.14 

Mb/Ma 

0.3877 

0.0060 


Table 3 

Instrument Parameters For Q Factor Calculation 


Instrument 

Resolution 

•^min 

(A) 

•^max 

(A) 

Eff. Aperture 
Radius (m) 

Total Throughput 
(%) 

Overhead 

(min) 

# Targets / Obs. 
(H ~ 13.5) 

HARPS-North^ 

115000 

3830 

6930 

3.58 

8 

3 

1 

Keck HIRES*’ 

55000 

5000 

6200 

10.0 

13 

3 

1 

SOPHIE'’ 

75000 

3820 

6930 

1.93 

4 

3 

1 

HEX'* 

30000 

4076 

7838 

9.2 

3 

3 

1 

APOGEE'’ 

22500 

15100 

17000 

2.12 

16 

15 

89 


^ Instrument parameters taken from http://www.tng.iac.es/instruiiients/harps/,, 

^ Instrument specs as reported in IJohnson et al. I 1120111) . Efficiency is taken to be 18% from http://www2.keck.hawaii.edu/ 
inst/hires/throughput .pdf, minus an additional 30% loss due to absorption from the iodine cell. 

^ Instrument parameters taken from http://www.obs-hp.fr/guide/sophie/sophie-info.html,, 

^ Instrument parameters are for the 316g5936 cross-disperser in R — 30000 mode using ThAr as a wavelength calibration. 
Although most planet work has been done using the R — 60000 mode, this resolution is not required to achieve 100 m s~^ 
RV precision. Efficiency taken from the HRS exposure time calculator http://het.as.utexas.edu/HET/hetweb/Instruments/ 
HRS/exp/exp_calc .html,, and is calculated at the center of the telescope’s observability track. 

® Instrument parameters taken from [Wilson et al. I ll2012i) . Effective aperture radius includes a 30% loss due to obstruction 
of the 2.5m diameter jGunn et al. ~i [2ooa . 
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Table 4 

TRILEGAL Parameters Used In KIC 3861595 Background EB Blend Probability Calculation 


Parameter 

Value 

Distance modulus resolution of Galaxy components 
IMF for single stars 

Binary fraction 

Binary mass ratios 

Extinction model 

Extinction model hz^dust 

Extinction model 

Extinction calibration at infinity 

Icr extinction dispersion 

Solar Galactocentric radius Rq 

Solar height above the disk zq 

Thin disk model 

Thin disk zq 

Thin disk to 

Thin disk a 

Thin disk hji d 

Thin disk radial cutoffs 

Thin disk (©) 

Thin disk SFR+AMR 

Thick disk model 

Thick disk hz^td 

Thick disk hii,td 

Thick disk radial cutoffs 

Thick disk 0.^^ (©) 

Thick disk SFR+AMR 

Halo model 

Halo 

Halo qh 

Halo (©) 

Halo SFR+AMR 

Bulge model 

Bulge am 

Bulge ao 

Bulge y/x axial ratio 77 

Bulge zjx axial ratio ^ 

Bulge Sun-GC-bar angle 

Bulge (GC) 

Bulge SFR+AMR 

0.1 mag 

Chabrier lognormal 

0.3 

0.7 to 0.1 

Exponential disk of form exp (— \z\ /h^ dust) * 6 xp {—R/hji dust) 

110 pc 

100000 pc 

0.0378 

0 . 

8700 pc 

24.2 pc 

squared hyperbolic secant 

94.6902 pc 

5.55079E9 yr 

1.6666 

2913.36 pc 

0, 15000 pc 

55.4082 Mq pc-^ 

2-step SFR + Fuhrman’s AMR + a enhancement with age(yr) = 0.7350974 -I- 0 
squared hyperbolic secant 

800 pc 

2394.07 pc 

0, 15000 pc 

0.001 Mq pc~^ 

11- 12 Gyr const. SFR + Z=0.008 with cr[M/H] = 0.1 dex with age(yr) = 4 -|- 0 

Oblate spheroid 

2698.93 pc 

0.583063 

0.000100397 Mq pc-3 

12- 13 Gyr -|- Ryan & Norris [M/H] distribution with age (yr) = 4 -|- 0 
triaxial bulge 

2500 pc 

95 pc 

0.68 

0.31 

15° 

406 Mq pc~^ 

10 Gyr, Zoccali et al. 2003 [M/H] + 0.3 dex with age(yr) = 4 — 2E9 






